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Abstract 

We present a high precision lattice calculation of the average up/down, strange 
and charm quark masses performed with Nf = 2 twisted mass Wilson fermions. 
The analysis includes data at four values of the lattice spacing and pion masses as 
low as ~ 270 MeV, allowing for accurate continuum limit and chiral extrapolation. 
The strange and charm masses are extracted by using several methods, based on 
different observables: the kaon and the r] s meson for the strange quark and the D, 
D s and r] c mesons for the charm. The quark mass renormalization is carried out 
non-perturbatively using the RI-MOM method. The results for the quark masses 
in the MS scheme read: fn ud (2 GeV) = 3.6(2) MeV, m s (2 GeV) = 95(6) MeV and 
m c {rfi c ) = 1.28(4) GeV. We also obtain the ratios m s /m U( i = 27.3(9) and m c /m s = 
12.0(3). 



1 Introduction 



A precise knowledge of the values of the quark masses is of great importance for testing 
the Standard Model (SM) of particle physics. From a phenomenological point of view, 
several useful observables to constrain the SM or to search for New Physics depend on 
quark masses, thus requiring accurate quark mass values in order to allow for significant 
theory/experiment comparisons. From a more theoretical side, explaining the quark mass 
hierarchy, which is not predicted by the SM, is a deep issue and a great challenge. Lattice 
QCD calculations play a primary role in the determination of quark masses. Recently, the 
progress achieved thanks to several high statistics unquenched simulations is leading to a 
significant reduction of the uncertainty on the quark mass values pQ- [6]. 

In this letter we present an accurate determination of the average up/down, strange 
and charm quark masses performed by the European Twisted Mass Collaboration (ETMC) 
with Nf = 2 maximally twisted mass Wilson fermions. The high precision of this analysis 
is mainly due to the extrapolation of the lattice results to the continuum limit, based on 
data at four values of the lattice spacing, to the well controlled chiral extrapolation, which 
uses simulated pion masses down to M n ~ 270 MeV, and to the use of the non-perturbative 
renormalization constants calculated in [7j. The only systematic uncertainty which is not 
accounted for by our results is the one stemming from the missing strange and charm quark 
vacuum polarization effects. Those are not accessible to us with Nf = 2 flavor simulations. 
However, a comparison of Nf = 2 results for the up/down and strange quark masses to 
already existing results from Nf = 2 + 1 quark flavor simulations [8] indicates that, for 
these observables, the error due to the partial quenching of the strange quark is smaller 
at present than other systematic uncertainties. The same conclusion is expected to be 
valid for the effects of the strange and charm partial quenching in the determination of the 
charm quark mass. 

In this work, the calculation of the isospin averaged up/down quark mass, based on the 
study of the pion mass and decay constant, closely follows the strategy of [9]. At variance 
with the latter, however, in the present analysis we include data at four values of the lattice 
spacing, and use the same lattice setup for all quark mass analyses. 

For the strange quark mass, the main improvement with respect to our previous 
work [10J, which used data at a single lattice spacing only, is the continuum limit. More- 
over, the chiral extrapolation is performed by using either SU(2)- or SU(3)-Chiral Pertur- 
bation Theory (ChPT). In order to extract the strange quark mass we have used both the 
kaon mass and the mass of the (unphysical) t] s meson composed of two degenerate valence 
strange quarks. In both cases, the ultimate physical input is the kaon mass, together with 
the pion mass and decay constant. 

For the charm quark mass, similarly to the strange quark, we have investigated several 
experimental inputs to extract its value: the mass of the D, D s and r] c mesons. In the 
charm quark sector discretization effects require some care and having data at four lattice 
spacings helps in performing the continuum limit. 

The results that we obtain for the quark masses are, in the MS scheme, 

m ud (2 GeV) = 3.6(1)(2) MeV = 3.6(2) MeV, 
m s (2 GeV) = 95(2)(6) MeV = 95(6) MeV, 
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m c (m c ) = 1.28(3)(3) GeV = 1.28(4) GeV, 



(1) 



where the two separate errors are, respectively, statistical and systematic. We also obtain 
for the ratios of quark masses the values 



The calculation is based on the Nf = 2 gauge field configurations generated by the ETMC 
with the tree-level improved Symanzik gauge action [TT] and the twisted mass quark ac- 
tion [12] at maximal twist, discussed in detail in [9], [13] -[16]. We simulated Nf = 2 dy- 
namical quarks, taken to be degenerate in mass, whose masses are eventually extrapolated 
to the physical isospin averaged mass of the up and down quarks. As already mentioned, 
the strange and charm quarks are quenched in the present calculation. 

The use of the twisted mass fermions turns out to be beneficial, since the pseudoscalar 
meson masses, which represent the basic ingredient of the calculation, are automatically 
improved at 0(a) [17]. As discussed in [IH1 [HI US] , we implement non-degenerate valence 
quarks in the twisted mass formulation by formally introducing a twisted doublet for each 
non-degenerate quark flavor. In the present analysis we thus include in the valence sector 
three twisted doublets, (u, d), (s,s f ) and (c, c'), with masses /i s and /i c , respectively. 
Within each doublet, the two valence quarks are regularized in the physical basis with 
Wilson parameters of opposite values, r = — r' = 1. Moreover, we only consider in the 
present study pseudoscalar mesons composed of valence quarks regularized with opposite 
r. This choice guarantees that the squared meson mass m 2 PS differs from its continuum 
counterpart only by terms of 0(a 2 fi q ) and 0(a 4 ) [20l f2Tj . 

Details of the ensembles of gauge configurations used in the present analysis and the 
values of the simulated valence quark masses are collected in Tables [1] and [21 respectively. 
In order to investigate the properties of the various light, strange and charmed mesons, we 
simulate the sea and valence up/down quark mass in the range 0.15 vnP^ ys < fii < 0.5 m phys , 
the valence strange quark mass within 0.8 < fi s < 1.5™^, and the valence charm 

quark mass within 0.9 ^ fi c ^ 2.0mf ws , with m p s hys and m phys being the physical 

strange and charm masses. Quark propagators with different valence masses are obtained 
using the so called multiple mass solver method J22j [23] , which allows to invert the Dirac 
operator for several valence masses at a relatively low computational cost. 

The statistical accuracy of the meson correlators is significantly improved by using 
the so-called "one-end" stochastic method, implemented in [23], which includes all spatial 
sources at a single timeslice. Statistical errors on the meson masses are evaluated using 
the jackknife procedure. With 16 jackknife bins for each configuration ensemble we have 
verified that autocorrelations are well under control. Statistical errors on the fit results 
which are based on data obtained from independent ensembles of gauge configurations are 



m s /m ud = 27.3(5) (7) = 27.3(9) 
m c /m s = 12.0(3) , 



(2) 



which are independent of both the renormalization scheme and scale. 



2 Simulation Details 
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Table 1: Details of the ensembles of gauge configurations used in the present study: value of 
the gauge coupling (3; value of the lattice spacing a; lattice size V = L 3 xT in lattice units; 
bare sea quark mass in lattice units; approximate value of the pion mass; approximate 
value of the product m n L; number of independent configurations N c f 9 . 
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Table 2: Values of simulated bare quark masses in lattice units for each configuration 
ensemble in the light, strange and charm sectors. 

evaluated using a bootstrap procedure, with 100 bootstrap samples, which properly takes 
into account cross-correlations. 

The analysis is based on a study of the dependence of meson masses on renormalized 
quark masses, with data at the four simulated values of the lattice spacing simultaneously 
analyzed. For the quark mass renormalization constants = Zp 1 we use the results 
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obtained in [TJ, which read 



Z P \ = {0.411(12), 0.437(7), 0.477(6)} at = {3.8, 3.9, 4.05} , (3) 

in the MS scheme at 2 GeV. The errors given in eq. ([3]) are those quoted in [7J and 
do not account neither of discretization errors nor for the uncertainty associated to the 
perturbative conversion from the RI-MOM to the MS scheme. The former are taken care 
by performing on the renormalized quark masses the extrapolation to the continuum limit. 
The uncertainty associated to the conversion from the RI-MOM to the MS scheme is 
included in our final estimate of the systematic error on the quark masses and will be 
discussed in sect. 3. For the renormalization constant at j3 — 4.2, not calculated in [7], we 
use the preliminary result Zp S (2 GeV)| 4 . 2 = 0.501(20), where the conservative uncertainty 
is due to the preliminarity of the result. 

The uncertainty on Zp has been taken into account by including in the definition of 
the x 2 to be minimized in the fits a term 



(z P (a) - Z P (a))' 
5Zp(a) 2 



(4) 



for each value of the lattice spacing and for each bootstrap sample, where Z P [a) ± 5Zp(a) 
is the input value for the renormalization constant at the lattice spacing a and for the 
bootstrap i, and Z p {a) the corresponding fit parameter. This procedure corresponds to 
assuming for the renormalization constant a Bayesian gaussian prior [91 [25]. 

The simultaneous analysis of data at different values of the lattice spacing also requires 
the data conversion from lattice units to a common scale. For the analysis in the pion 
sector, we have expressed all dimensionful quantities in units of the Sommer parameter 
To [26]. We use for r^/a in the chiral limit the values 



ro 
a 



= {4.54(7), 5.35(4), 6.71(4), 8.36(6)} at = {3.8, 3.9, 4.05, 4.2} , (5) 

/3 



obtained from an extension of the analyses in 0, EH] with the inclusion of all four lattice 
spacings. As in [9], the chiral extrapolation of r^/a is performed by using three ansatze 
for the sea quark mass dependence: linear only, quadratic only and quadratic+linear. The 
size of mass- dependent discretization effects is verified by including in the fits 0{a 2 m{) and 
0(a 2 mf) terms, which turn out to be negligible. The uncertainties on the results given in 
eq. flSJ) include the systematic errors estimated as the spread among the values obtained 
from the above-mentioned fits. In the present analysis the uncertainty on the r /a values 
is taken into account by adding a term to the \ 2 of the fit in a similar way to Zp, as 
explained above. 

The analysis in the pion sector is also used to determine, besides the value of the average 
up/down quark mass at the physical point, the lattice spacing at each coupling (3. The 
physical input used for this determination is the pion decay constant f w . In the successive 
determination of the strange and charm quark masses, the data are analyzed directly in 
physical units. 
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3 Up/down quark mass 



The calculation of the up/down quark mass follows the strategy of [9J. The analysis is 
repeated here including simultaneously all data available at the four values of the lattice 
spacing. 

We have studied the dependence of the pion mass and decay constant on the renor- 
malized quark mass. For these quantities the predictions based on NLO ChPT and the 
Symanzik expansion up to 0(a 2 ) can be written in the form 



mi = (2 Bonn) 



2 B Q mi , f2B mi\ „ 9 ( „ „ , / 2 B Q mi 



2B m, (2B m l \ 2 / . . . 

I6^7f log (l^7fJ + 4m ' + a A 5 6 log|T ^ 11 ' (6) 




where mi is the renormalized light quark mass and Bq and fo are the low energy constants 
(LECs) entering the LO chiral LagrangiaijE 

The coefficients of the discretization terms of 0(a 2 mi log(m^)) receive a contribution 
from the 0(a 2 ) splitting between the neutral and charged pion (squared) mass, Am 2 = 
m^o — ??7^± , which occurs with twisted mass fermions. This contribution has been recently 
evaluated in [27J. Our main results are obtained through a fit of eqs. (jSJ), with the coeffi- 
cients P3 and Pq obtained by expanding the results of [27] up to 0(a 2 ). We have also verified 
that the results obtained in this way are indistinguishible from those obtained using directly 
the resummed formulae of |27j. From our fit, the splitting Am 2 turns out to be deter- 
mined with an uncertainty of approximately 60%. We obtain Am 2 = — (33 ± 19) a 2 Aq CD , 
which is consistent with a direct ETMC determination performed with two lattice spacings 
(Am 2 = —50 a 2 Aq CD [28]). On the final result for the light quark mass the impact of this 
correction is at the level of the fitting error. 

Lattice results for pion masses and decay constants have been corrected for finite size 
effects (FSE) evaluated using the resummed Liischer formulae. The effect of the 0(a 2 ) 
isospin breaking has been taken into account also in these corrections [29]. On our pion 
data, FSE vary between 0.2% and 2%, depending on the simulated mass and volume. The 
inclusion of the pion mass splitting in the FSE induces a variation of about 15-40% in the 
finite size correction itself. This effect is at the level of one third of the statistical error 
for our lightest pion mass at /3 = 3.9 on the smaller volume, and even smaller in the other 
cases. 

In fig. [1] (left) we show the dependence of r$ m 2 /mi on the renormalized light quark mass 
at the four ft's, and the curves corresponding to the best fit of the lattice data according 
to eq. (|6]). 

In order to illustrate the dependence of the pion mass on the lattice cutoff, we have 
interpolated the lattice data for m 2 at the four values of the lattice spacing to a common 
reference value of the light quark mass, rrf^J = 50 MeV. The resulting values of (r m n ) 2 
obtained in this way are shown in fig. [1] (right) as a function of a 2 , together with the 



1 The pseudoscalar decay constant in the chiral limit, fo, is normalized such that f„ — 130.7 MeV at 
the physical point. 
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Fig. 1: Left: Dependence of r^m^/mi on the renormalized light quark mass at the four 
lattice spacings. Right: Dependence of (r Q m n ) 2 on the squared lattice spacing, for a fixed 
reference light quark mass (n\J = 50 MeVj. Empty diamonds represent continuum limit 
results. 

corresponding continuum extrapolation. We see that discretization errors on the pion 
mass square are below 10% at = 3.9 and negligible within the fitting errors at /3 = 4.2. 

The value of the physical up/down quark mass is extracted from the ratio m 2 /f 2 
using as an input the experimental value of the latter ratio|§ In order to estimate the 
systematic uncertainty due to discretization effects we have performed both a fit without 
the logarithmic discretization terms, i.e. with P 3 = P 6 = in eq. ([6]) (the so called fit 
B of [9]), and a fit without all 0(a 2 ) corrections, i.e. with P2 = P3 = P5 = Pq = 
(the so called fit A of [9]). Both these ansatze turn out be compatible with the lattice 
data. We find that the result for the up/down quark mass decreases by approximately 2% 
and increases of about 6% in the two cases respectively, so that we estimate an overall 
uncertainty due to residual discretization effects of ±4%. We have also tried to add in the 
fit discretization terms of 0(a 2 m 2 ) or (9(a 4 ). In both cases these terms turn out to be 
hardly determined with our data, leading for m u d to results consistent with those obtained 
from the other fits, but with uncertainties larger by a factor three. 

For estimating the systematic uncertainty due to the chiral extrapolation we have also 
considered a fit including a NNLO local contribution proportional to the light quark mass 
square. In this case we are not able to determine all the fitting parameters and we are thus 

2 In order to account for the electromagnetic isospin breaking effects which are not introduced in the 
lattice simulation, we have used in the present analysis as "experimental" values of the pion and kaon mass 
the combinations [30] 

{Ml) QCD = Mlo , {M 2 K ) QCD = \ [M 2 K0 + M| + - (f + A E )(M^ + - M^)] , (7) 

with = f. The values of the experimental inputs for the pion and kaon masses are then m% xp = 
135.0 MeV, m exp = 494.4 MeV. 
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forced to introduce, on the additional LECs, priors as in [9]. In this way we find that the 
result for m u d increases by 6%. 

The results of the fits described above are collected in Table [3] in the appendix. 

As anticipated in the previous section, we also include in the final result a systematic 
uncertainty coming from the perturbative conversion of the quark mass renormalization 
constant from the RI-MOM to the MS scheme. Using the results of the 3-loop calculation 
of [31], one can write the relation between the quark mass in the two schemes as 

; " (/,) - 1 - 0.424a s (/i) - 0.827a s (/i) 2 - 2.126a s (/i) 3 + 0(a s (fi) 4 ) . (8) 



m RI (/i) 

The uncertainty due to the truncation of the perturbative series has been conservatively 
estimated by assuming the unknown 0(a^) term to be as large as the 0(a^) one. Evaluating 
this term at the renormalization scale /i ~ 3 GeV, which is the typical scale of the non- 
perturbative RI-MOM calculation in our simulation [TJ, and using a s (3 GeV,Nf = 2) = 
0.202, we then find that this uncertainty corresponds to approximately ±2 %. 

Adding in quadrature the three systematic errors discussed above we obtain m u d = 
3.55(14) (ill) MeV in the MS scheme at the renormalization scale of 2 GeV, where the 
two errors are statistical and systematic, respectively. Finally we symmetrize the error, 
obtaining 

m ud {2 GeV) = 3.6(1) (2) MeV = 3.6(2) MeV. (9) 

Note that, in the symmetrized result, the uncertainties due to discretization effects, chiral 
extrapolation and perturbative conversion give similar contributions to the final systematic 
error, at the level of 4%, 3% and 2% respectively. 

Using as an input the experimental value of the pion decay constant, the fits also provide 
us with the values of lattice spacing at the four simulated /3's, which are used in the rest 
of the analysis. They read, at (3 = {3.8, 3.9,4.05,4.2} respectively, 

a\f) = {0.098(3)(2), 0.085(2)(1), 0.067(2)(1), 0.054(1)(1)} fm , (10) 

where again the two errors are statistical and systematic. The results in eqs. © and (fit)]) 
are in good agreement with the previous ETMC determination obtained in [9] from the 
analysis of data at (3 = 3.8, 3.9 and 4.05. 

We observe that, in principle, the ratio of lattice spacings at two different values could 
be determined from the fit of the pion meson mass and decay constant, without using the 
additional information coming from the values of r^/a of eq. (jSJ). With our data, however, 
the uncertainties on the values of the quark mass renormalization constant, as well as the a 
priori unknown size of discretization errors affecting the pion masses and decay constants, 
do not allow to achieve a reliable determination of these ratios. 



4 Strange quark mass 

In this section, we first present the determination of the strange quark mass based on the 
study of the kaon meson mass. The alternative determination based on the study of the 
7] s meson will be discussed later on. 
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Since the valence strange quark mass has not been previously tuned in our simulation, 
the determination of the physical strange quark mass requires an interpolation of the lattice 
data. As already mentioned, for all values of the lattice couplings, the simulated values of 
the strange quark masses are approximately in the range 0.8™^ ^ ji s ^ l.bm p s hys . 

In order to better discriminate the strange quark mass dependence of the kaon masses 
on other dependencies, in particular discretization effects, we firstly slightly interpolate 
the lattice data to three reference values of the strange quark mass, which are chosen to 
be equal at the four lattice spacings: rrf/f = {80 , 95 , 110} MeV. The interpolations to 
the reference masses are performed by using quadratic spline fits. Then, at fixed reference 
strange mass, we simultaneously study the kaon mass dependence on the up/down quark 
mass and on discretization effects, thus performing the chiral extrapolation and taking 
the continuum limit. In this step, we have considered chiral fits based either on SU(2)- 
ChPT [31 [32], where the chiral symmetry is assumed for the up/down quark only, or 
partially quenched SU(3)-ChPT [33J, where instead also the valence strange quark is treated 
as light. In order to extrapolate the kaon mass values to the continuum and to the physical 
m u d limit, we use the results for the average up/down quark mass and for the lattice 
spacings obtained in eqs. (jHJ) and (fTD"]) . at each reference value of the strange quark mass. 
Finally, we study the kaon mass dependence on the strange quark mass, and determine the 
value of the physical strange quark mass using the experimental value of mjf. 

Let us describe the chiral fits in more detail. As discussed above, fits are performed 
in two steps: 1) the strange quark mass is fixed to the reference values and only the mi 
and a 2 dependence of the kaon mass is studied; 2) the so obtained data in the continuum 
limit and at the physical m u d value are studied as a function of the strange quark mass. 
In these two steps, we have considered for the kaon meson mass functional forms based on 
the predictions of either NLO SU(2)-ChPT [3J, which predicts the absence at this order of 
chiral logs, 

1) m 2 K (m s ,mi,a) = Q^ms) + Q 2 (m s )m l + Q 3 (m s ) a 2 , Vm s , (11) 

2) m 2 K (m s ,m 1 i hys ,a = 0) = Qi(m s ) + Q 2 (^ s ) mf 1 ^ = Q 4 + Q 5 m s , (12) 
or SU(3)-ChPT [33] 

1) m 2 K (m s ,mi,a) = B ■ (m s + mi) ■ (1 + Q e (m s ) + Q 7 {m s ) m t + Q%{m 8 ) a 2 ) , Vm s , 

(13) 

2) m 2 K (m s , mf hys , a = 0) = B ■ (m s + mf ys ) ■ (1 + Q 6 (m s ) + Q 7 (m s ) mf ys ) = 

„ / nhnQ^ I 2 Bn Trio 2 Bom* „ \ , . 

= B - (m s + mf ys )- U + TT^tJ log — ^-f + Q 9 m s , 14 

V (4vr/ ) 2 (4vr/ ) 2 / 

where B and /o are determined from the pion fit described in the previous section. Note 
that the dependence of the kaon mass on the strange quark mass is not determined by 
the chiral symmetry in the SU(2) theory. We find that, with our choice of three reference 
strange masses around the physical value, a linear fit as given in eq. ffl2|) is perfectly 
adequate to describe the data. 
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Fig. 2: Left: Dependence of m 2 K and m 2 s on the renormalized light quark mass, for a Gxed 
reference strange quark mass (rrf/* = 95 MeV ) and at the four lattice spacings. The orange 
vertical line corresponds to the physical up/down mass. Right: Dependence ofm 2 K and m 2 e 
on the squared lattice spacing, for a hxed reference strange quark mass (rrf s e f = 95 MeVj 
and at the physical up/down mass. Empty diamonds represent continuum limit results. 

For illustration we show in fig. [2] the combined chiral/continuum fit based on SU(2)- 
ChPT, eq. (TTTj) . for a fixed reference value of the strange quark mass, as a function of the 
light quark mass (left) and of the squared lattice spacing (right). In fig. [3] the dependence 
on the strange quark mass is shown, in the case of the SU(2) analysis (see eq. (fl2"]l ). The 
dependencies are shown for the kaon squared mass as well as for the n s squared mass 
discussed hereafter. 

As an alternative way to determine the strange quark mass we have studied the depen- 
dence on m s of a meson made up of two strange valence quarks j5]. The advantage of this 
approach is that the mass of this unphysical meson, denoted as rj s , is only sensitive to the 
up/down quark mass through sea quark effects, and it is thus expected to require only a 
very smooth chiral extrapolation. This expectation will be confirmed by our analysis. The 
price to pay is the need for an additional chiral fit required to determine the r] s mass at 
the physical point. 

In the real world, the n s meson is known to mix with the (uu + dd) component to 
produce the physical r\ and r/ mesons. This mixing proceeds through the contribution of 
disconnected diagrams, which are known to be rather noisy on the lattice and therefore 
computationally expensive. In order to avoid this computation we consider here the two 
strange quarks composing the meson as degenerate in mass but distinct in flavor. Though 
this n s meson does not exist in nature, its mass can be determined on the lattice [5]. 

In order to relate the mass of the n s meson to the physically observable and rrix, 
we have studied its dependence on the kaon and pion masses for different values of the 
simulated light and strange quark masses. This dependence turns out to be well described 
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m (2 GeV) [MeV] 



Fig. 3: Dependence of m\ and m 2 , in the continuum limit and at the physical up/down 
mass, on the strange quark mass. The strange mass results obtained from the SU(2)-ChPT 
analyses of kaon and rj s mesons are also shown, with empty diamonds. 



by both the functional fornix based on either NLO SU(2)-ChPT, 



'Is 

or SU(3)-ChPT 



mi =R 1 + R 2 (2m 2 K -ml) + R 3 ml + R 4 a 2 , (15) 



ml = (2 m 2 K -ml) - 1 + - £,) log(2 £,) + (R 7 + 1) - &) + R 8 a 

,2 



-ml [-£ log(2 &) + 6 log(2 6) + #7 (6 - 6)] > (16) 

with = ml/(A-n fo) 2 and £ s = (2m^ — ml)/(4n f ) 2 . We observe that, within the 
accuracy of our lattice data, the 0(a 2 ) term in the T] s mass is found to be independent of 
the strange quark mass. 

Once the physical values of the kaon and pion mass are inserted in eqs. (j!5p and (1161) . 
we find that the two fits yield very close results for the r) s meson mass, namely 

m Va = 692(1) MeV from SU(2) , m Va = 689(2) MeV from SU(3) , (17) 

to be compared with the LO SU(3) prediction (m Va ) LO = J 2 m\ — m\ = 686 MeV and 
with the lattice determination of |5] m Vs = 686(4) MeV. 

Once the mass of the n s meson has been determined, the strange quark mass can be 
extracted by following the very same procedure described for the case of the kaon mass. 
At first, lattice data at fixed reference strange mass are extrapolated to the continuum 
and to the physical up/down mass (see fig. |2]). After this extrapolation, the value of the 



3 The functional forms in eqs. (| 1 5 [) - (|TH]) are obtained from the ChPT formulae given in eqs. (IT51) - (j!?Tj) 
by replacing quark masses in terms of meson masses, and keeping terms up to NLO. 
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physical strange quark mass is extracted by studying the dependence on the strange mass 
(see fig. |3]). We have considered the following fitting functions based on NLO-ChPT for 
the dependence of the rj s meson 1) on the (sea) up/down quark mass and on the leading 
discretization effects, and 2) on the strange quark mass: 

1) m 2 g (m s ,m h a) =Ti(m s ) + T 2 (m s )m l + T 3 (m s )a 2 , Vm s , (18) 

2) m 2 rjs (m s ,m 1 f iy \a = Q) =Ti(m s ) + T 2 {m s ) mf 1 ^ = T 4 + T 5 m s , (19) 
in SU(2), and 

1) m 2 s (m 8 ,mi, a) = 2B m s ■ (1 + T 6 (m 8 ) + T 7 ( y m s )m l + T 8 (m s ) a 2 ) , Vm s , (20) 

2) m 2 g (m s , m^ hys , a = 0) = 2 B m s ■ (1 + T 6 (m s ) + T 7 (m s ) mf iys ) = 

= 2B m s - 2 log 2 +T 9 + T w m s , 21 

V ( 47r /o) 2 (4vr/ ) 2 / 

in SU(3). The LECs T 2 and T 7 , describing the dependence on the light quark mass, are 
found to be independent of the strange mass, within the accuracy of our lattice data. They 
are then fitted with a single parameter for all reference strange quark masses. 

The results for the strange quark mass obtained from both the kaon and the rj s meson 
masses turn out to be well consistent, as it can be seen from Table H]in the appendix where 
the several m s values obtained from different fits are collected. 

In order to quote a final estimate for the strange quark mass we choose as a central 
value the weighted average of the results from the four determinations discussed above, 
namely from K and 7] s and based on SU(2)- and SU(3)-ChPT. In the MS scheme at 2 GeV 
this average reads m s (2 GeV) = 95(2) MeV, with the 2 MeV error representing the typical 
statistical and fitting uncertainty. The difference between the determinations based on the 
K and 7] s mesons is about 3%. The results obtained from either the SU(2) or the SU(3) 
fits are practically the same in the analysis based on the rj s and differ by approximately 
3% in the kaon case. In order to evaluate the uncertainty of the continuum extrapolation 
we have proceeded in two ways. We have either added an 0(a 4 ) mass independent term 
in eqs. ( TTTj) . (fT3|) . ( ITS]) and (|20|) . or we have excluded from these fits the data from the 
coarser lattice (with a ~ 0.098 fm). We find that the 0(a 4 ) term turns out to be hardly 
determined in the fit, leading to a factor three larger uncertainties. The exclusion of (3 = 3.8 
data, instead, yields a variation of the results of approximately 2% leaving the fitting error 
approximately unchanged. We then assume ±2% as uncertainty related to the continuum 
extrapolation. The different fits considered for the determination of the up/down mass 
and of the lattice spacing affect the determination of the strange mass at the level of 3%. 
Finally, we include also in this case an uncertainty of 2% related to the truncation of the 
perturbative expansion in the conversion from the RI-MOM to the MS scheme. Combining 
all these uncertainties in quadrature, we quote as our final estimate of the strange quark 
mass in the MS scheme 

m s (2 GeV) = 95(2)(6) MeV = 95(6) MeV. (22) 
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We observe that our result for the strange mass in eq. (1221) is, though compatible, smaller 
than the value obtained in [10] at a fixed value of the lattice spacing (a ~ 0.085 fm). 
This is a consequence of discretization effects, which are at the level of 15% in m\ on 
the a ~ 0.085 fm lattice, as shown in fig. |2] (right). A further comparison can be done 
with the ETMC estimate of the strange quark mass that appeared in the recent work 
on the bag parameter [34]. Within that analysis, based on data at three ft values 
(3.8, 3.9 and 4.05), the strange quark mass is determined from the same lattice setup by 
performing an SU(2) chiral fit of the kaon meson mass. The result obtained in [34] reads 
m s (2 GeV) = 92(5) MeV, to be compared to our result m s (2 GeV) = 92.1(3.8) MeV, 
obtained from the same fit (see Table H]). 

Using our determinations of both the strange and light quark masses, we can obtain a 
prediction for the ratio m s /m u d, which is both a scheme and scale independent quantity. 
The several m s /m u d values obtained from different fits are collected in Table [5] in the 
appendix]^ The result that we quote as our final estimate is 

m s /m ud = 27.3(5) (7) = 27.3(9). (23) 

5 Charm quark mass 

The determination of the charm quark mass follows, quite closely, the strategy adopted 
in the determination of the strange quark mass discussed in the previous section. In this 
case, we use as experimental input the masses of the D, D s and rj c mesons. 

As for the strange quark case, the analysis requires an interpolation of the lattice data, 
being the simulated charm masses roughly in the range 0.9mf 9S ^ /i c ~ 2.0m^ hys . In 
order to better study the a 2 and mi dependence of charmed meson masses, we first use 
a quadratic spline fit to interpolate the data at three reference values of the charm mass 
which are equal at the four /3 values: frf c eJ \2 GeV) = {1.08, 1.16, 1.24} GeV. We have 
verified that a different choice of the values of the reference masses leaves the charm quark 
results unchanged. At fixed reference charm mass, we then study the dependence of the 
D, D s and r] c meson on the up/down mass (and on the strange mass in the case of the 
D s meson) and on discretization terms, thus getting the results for the meson masses 
in the continuum limit, at the physical values of the light (and strange) quark masses, 
and at the reference charm mass. Finally, the value of the physical charm quark mass is 
extracted by fitting these data as a function of the charm quark mass and using as an 
input the experimental value of the corresponding charmed meson m e ^ p = 1.870 GeV, 
rrfiP = 1.969 GeV, m£* = 2.981 GeV@ 

4 The results for the ratio m s /m u d collected in Table [5] are slightly different from the ratios of the 
m s and m u d results. This difference originates from the fact that in the ratio m s /m u d the quark mass 
renormalization constant Zp 1 exactly cancels out, whereas in the determinations of m s and ra u d the central 
values of Zp are slightly modified by the fitting procedure. 

5 The experimental value of the meson masses should be corrected to take into account the absence of 
electromagnetic effects and, in the case of the r) c , of disconnected diagrams in the lattice calculation. For 
the r\ c meson, these corrections are estimated to be of the order of 5 MeV [5] , thus affecting the extracted 
charm quark mass to approximately 0.2%. Similar corrections are expected for the D and D s mesons. 
Given our uncertainties, we can safely neglect these corrections in the analysis. 
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In order to fit the meson masses we have considered the following (phenomenological) 
polynomial fits, which turn out to describe well the dependence on the light and strange 
quark masses and on the lattice cutoff of the D, D s and rj c meson masses, at fixed (reference) 
charm mass m c , 

m H (m c ,m s ,mi,a) = Cf(m c ) + Cf (m c ) raj + Cf(ra c )ra s + Cf(ra c ) a 2 , Vra c , (24) 

with H = D t D s , r] c . From the fits, we find that the coefficients Cf and Cf turn out to be 
independent of the charm mass within the statistical errors. The latter coefficient Cf , of 
course, enters the fit only in the D s case. 

For the charm mass dependence, a constant plus either a linear or a l/m c term have 
been considered for describing data of the D, D s and r\ c mesons, namely 

m H (m c , m p s hy % mf y \ a = 0) = Cf (m c ) + Cf (m c ) m phys + Cf (m c ) mf ys = 

= Cf + Q- + C?m c . (25) 

m c 

Since we have data at three reference charm masses (close to the physical charm), we can 
keep only one of the coefficients Cf , Cf different from zero. We find that both choices 
describe very well the lattice data and affect only in a marginal way the interpolation to 
the physical charm mass. 

In fig. H]we show the dependence of the D, D s and r] c masses on the light quark mass 
at a fixed reference charm mass, for the four /3's. For the D s and r\ c mesons, which contain 
the light quark in the sea only, this dependence turns out to be practically invisible. 

In fig. [5] (left) the meson masses at physical light and strange quark masses are shown 
as a function of a 2 , for a reference value of the charm quark mass. As can be seen from this 
plot, discretization effects on the r\ c meson mass vary from approximately 4% on the finest 
lattice up to 14% on the coarsest one. These effects are larger than those affecting the D 
and D s meson masses by approximately 30%. Fig. [5] (left) also shows that the dependence 
of the three charmed meson masses on a 2 is very well described by a linear behaviour, and 
attempts to vary the continuum extrapolation with respect to the simple linear fit produce 
only small effects. The latter are included in the estimate of the systematic uncertainty, 
as discussed below. 

Finally, fig. [5] (right) shows the dependence of the D, D s and i] c masses on the charm 
mass (in the continuum limit and at physical light and strange mass) and the interpolation 
to the physical charm. 

Using the experimental values of the considered meson masses we get for the charm 
quark mass the results collected in Table El 

In order to estimate the uncertainty due to the continuum extrapolation we have pro- 
ceeded in two ways. We have either added in the fitting form of eq. fl24j) an 0(a 4 ) depen- 
dence, which turns out to be hardly determined thus leading to uncertainties larger by a 
factor three, or we have excluded the data from the coarser lattice (with a ~ 0.098 fm). 
This latter analysis yields a variation of the results of approximately 1.5%. The two de- 
pendencies of the meson masses on the charm quark mass, considered in eq. (!25j) . yield 
results that differ by only few MeV. The systematic uncertainty then comes from the sum 
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Fig. 4: Left: Dependence of (left) and mo s and m Vc (right) on the light quark mass, 
at fixed reference charm quark mass (fff c e f = 1.16 GeV) and for the four simulated lattice 
spacings. For the D s meson the strange quark mass is fixed to the reference value nf s e * = 
95 MeV. 

in quadrature of the approximately 1% spread among the three determinations from the 
D, D s and rj c mesons, the 1.5% uncertainty due to discretization effects and the 2% uncer- 
tainty coming from the the perturbative conversion of the renormalization constants from 
the RI-MOM to the MS scheme. 

We quote as our final result for the charm quark mass in the MS scheme 

m c (2GeV) = 1.14(3) (3) GeV = 1.14(4) GeV (26) 
-> m c (m c ) = 1.28(4) GeV, 

where the evolution to the more convential scale given by fn c itself has been performed at 
N 3 LO [31] with N f = 2, consistently with our non-perturbative evaluation of the renormal- 
ization constant. Had we evolved with Nf = 4, which is the number of active flavours above 
\i = m c , the result for m c (m c ) would have increased by less than one standard deviation. 

Our result is compatible with the preliminary estimate of the charm quark mass, 
m c (2 GeV) = 1.23(6), obtained by ETMC [35] using data at three lattice spacings and 
preliminary values for the renormalization constants. It is also in good agreement with 
the HPQCD result m c (m c ) = 1.268(9) GeV [6], with a larger uncertainty in our determi- 
nation. Finally, our result is in good agreement with the recent sum rules determination 
m c (m c ) = 1.279(13) GeV of [36]. 

We also provide a prediction for the scheme and scale independent ratio m c /m s . The 
several m c j m s values obtained from different fits are collected in Table [7] in the appendix. 
The result that we quote as our final estimate is 

m c /m s = 12.0(3) , (27) 

in good agreement with the other recent lattice determination m c /m s = 11.85(16) [5]. 
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Fig. 5: Left: Dependence of mo, rnD a and m Vc , at fixed reference charm quark mass 
(jffzf — 1.16 GeV) and at physical up/down and strange quark mass, on the squared 
lattice spacing. Empty diamonds represent continuum limit results. Right: Dependence of 
rri£j, m£) s and m Vc , in the continuum limit and at physical up/down and strange quarks, 
on the charm quark mass. The charm mass results from the three determinations are also 
shown, with empty diamonds. 



6 Conclusions 

We have presented results for the average up/down, strange and charm quark masses, ob- 
tained with Nf = 2 twisted mass Wilson fermions. The analysis includes data at four values 
of the lattice spacing and pion masses as low as ~ 270 MeV, allowing a well controlled 
continuum limit and chiral extrapolation. Within the strange sector the chiral extrapola- 
tion is performed by using either SU(2)- or SU(3)-ChPT. The strange and charm masses 
are extracted by using several methods, based on different meson mass inputs: the kaon 
and the rj s meson for the strange quark and the D, D s and r\ c mesons for the charm. The 
quark mass renormalization is carried out non-perturbatively using the RI-MOM method. 

The results for the quark masses in the MS scheme read: rn u d{2 GeV) = 3.6(2) MeV, 
m s (2 GeV) = 95(6) MeV and m c (rn c ) = 1.28(4) GeV. The quoted errors include the 
uncertainty in the perturbative conversion of the renormalization constants from the RI- 
MOM to the MS scheme, which is conservatively estimated to be at the level of 2%. 
We emphasize that this uncertainty is not related to the lattice calculation itself, but 
comes from continuum perturbation theory. If the RI-MOM scheme was chosen as a 
reference scheme and, say, 3 GeV as a reference scale, which is the typical scale of the 
non-perturbative RI-MOM calculation in our lattice simulation [7] , this uncertainty would 
not be present at all. For reference we provide our results for the quark masses also in this 
scheme: m£[(3 GeV) = 3.9(1) (2) MeV, mf (3 GeV) = 102(2) (6) MeV and m* : (3 GeV) = 
1.22(3)(2) GeV. 

We have also evaluated the quark mass ratios m s /m u d = 27.3(9) and m c /m s = 12.0(3), 
which are independent on both the renormalization scale and scheme. 

The only systematic uncertainty which is not accounted for by our results is the one 
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stemming from the missing strange and charm quark vacuum polarization effects. A com- 
parison, for instance, of our Nf = 2 result for the strange quark mass, to already existing 
results from Nf = 2 + 1 quark flavor simulations [8] indicates that the error due to the 
partial quenching of the strange quark is smaller at present than other systematic uncer- 
tainties. In this respect we mention that simulations with Nf = 2 + 1 + 1 dynamical flavors 
are already being performed by ETMC and preliminary results for several flavor physics 
observables have been recently presented [371 EH] • 

We thank all the ETMC members for fruitful discussions and the apeNEXT computer 
centres in Rome for their invaluable technical help. Some computation time has been 
used for that project on the BlueGene system at IDRIS. We are also grateful to Gilberto 
Colangelo for having provided us with a routine for the calculation of the FSE [29], and to 
Chris Sachrajda for drawing our attention on the uncertainty in the perturbative conversion 
of the renormalization constants. 

7 Appendix 

We collect in this appendix the results obtained for the up/down, strange and charm quark 
masses from the different fits considered in the present analysis. 

As discussed in sect. 3, in the light quark sector we have performed the following fits: 

- LI: this is our best fit which is based on NLO ChPT with the inclusion of 0(a 2 ) 
discretization effects. This fit corresponds to eq. (jSJ) with all parameters different 
from zero. 

- L2: same as LI but without discretization terms, i.e. P2 = P3 = P5 = Pq = in 
eq. ©. 

- L3: same as LI with the inclusion of a NNLO correction proportional to the square 
of the light quark mass. 

The results for the up/down quark mass obtained from these fits are collected in Table [31 
For illustration, we also show in the table the value of m u d obtained from a fit (denoted as 
L4 here and B in [9]) without logarithmic discretization terms, i.e. with P3 = P$ = in 
eq. (jSJ), and without isospin breaking corrections in the FSE. 





LI 


L2 


L3 


L4 


mi [ MeV] 


3.55(14) 


3.75(7) 


3.78(17) 


3.47(11) 



Table 3: Results for the up/down quark mass in the MS scheme at 2 GeV, as obtained 
from fits LI, L2, L3 and L4. 

For the strange quark mass we collect the results of the different fits in Table HI where 
we use the short notation: K-SU(2), K-SU(3), ?7 S -SU(2), and r] s -SU(3) for distinguishing 
the determinations from the kaon and rj s masses, and based on SU(2)- and SU(3)-ChPT. 
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m s [ MeV] 


K-SU(2) 


K-SU(3) 


r/,-SU(2) 


%-SU(3) 


LI 


92.1(3.8) 


94.7(2.2) 


96.0(2.6) 


95.5(2.1) 


L2 


91.6(3.9) 


94.6(2.3) 


95.4(2.6) 


95.3(1.9) 


L3 


95.4(3.8) 


94.7(2.1) 


99.4(2.9) 


97.7(2.2) 



Table 4: Results for the strange quark mass in the MS scheme at 2 GeV, as obtained from 
the different fits within the light and strange quark sectors. 



m s /m ud 


K-SU(2) 


K-SU(3) 


»fc-SU(2) 


r? s -SU(3) 


LI 


26.9(5) 


27.2(5) 


27.6(4) 


27.3(7) 


L2 


27.1(5) 


26.9(3) 


27.5(3) 


26.8(3) 


L3 


25.7(5) 


26.0(6) 


26.5(6) 


26.0(7) 



Table 5: Results for the ratio m s /m U( i, as obtained from the different fits within the light 
and strange quark sectors. 

The results for the ratio m s /m u d are given in Table El 

Finally, for the charm quark mass and the ratio m c /m s the results are collected in 
Tables HD and [7] , where D, D s or r] c indicate the meson whose mass is used as input. These 
analyses are practically insensitive to the choice of the fit in the pion sector and only the 
results obtained from the fit LI are shown in the tables. Similarly, for the ratio m c /m s the 
values shown in Table [7] correspond to the analysis of the D meson only, since the analyses 
of the D s or rj c mesons yield practically identical results. 



D 


D s 




m c [GeV] 1.14(3) 


1.14(3) 


1.15(2) 



Table 6: Results for the charm quark mass in the MS scheme at 2 GeV, as obtained from 
the different fits within the charm sector. 





K-SU(2) 


K-SU(3) 


^ S -SU(2) 


r/ s -SU(3) 


m c /m s 


12.4(4) 


12.1(2) 


11.9(2) 


12.0(3) 



Table 7: Results for the ratio m c /m s , as obtained from the different fits within the strange 
quark sector, and from the analysis of the D meson mass in the charm sector. 
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